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Abstract 

The results of the MACRO experiment on atmospheric neutrino oscillations are 
presented and discussed. The data concern different event topologies with average 
neutrino energies of ~ 4 and ~ 50 GeV. The Multiple Coulomb Scattering of muons 
in the MACRO absorbers was used to estimate the neutrino energy of each event 
of the higher energy sample. The angular distributions, the absolute fluxes and 
the L/E distribution strongly favour Vn — > v T oscillations with maximal mixing 
and Am 2 = 0.0025 eV 2 . Results are presented on the searches for astrophysical 
sources of high energy muon neutrinos, for bursts of electron antineutrinos from 
stellar gravitational collapses and on indirect searches for WIMPs from the Earth 
and from the Sun. 



1 Introduction 



MACRO was a large area multipurpose underground experiment designed to search for 
rare events in the penetrating cosmic radiation. These included the study of atmospheric 
neutrinos and their oscillations, high energy (E u £ 1 GeV) muon neutrino astronomy, 
indirect searches for WIMPs and search for low energy (E u ^ 7 MeV) stellar collapse 
neutrinos. The detector was placed in the Gran Sasso laboratory, located on the highway 
Rome-Teramo , 120 km east of Rome. The lab consists of three underground halls, each 
about 100 m long. It is at an altitude of 963 m above sea level, is well shielded from 
cosmic rays by a mean rock thickness of ~ 3700 m.w.e.; the minimum is 3150 m.w.e: this 
defines the minimum muon energy at the surface at ~ 1.3 TeV in order to reach MACRO. 
The average residual muon energy and the muon flux are ~ 320 GeV and ~ lm _2 h _1 , 
respectively. 

1 MACRO Collaboration. See Ref. 1 for a list of MACRO Collaborators and Institutions. 
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Figure 1 : Cross section of the MACRO detector and sketch of different event topologies induced 
by Vp, — > A* charged current interactions. The black points and the black rectangles represent 
streamer tubes and scintillator hits, respectively. Tracking was performed by the streamer tubes; 
the time-of-flight (T.O.F.) of the muons was measured by the scintillators for Up Semicontained 
(Internal upgoing - IU /i) and Upthroughgoing events (and also for downgoing muons). 



MACRO was composed of three sub-detectors: liquid scintillation counters, limited 
streamer tubes and nuclear track detectors. Each one of them could be used in "stand- 
alone" and in "combined" mode. A cross section of the detector is shown in Fig. 1. Notice 
the division in the lower and in the upper part (this was often referred to as the Attico); 
the inner part of the Attico was empty and lodged the electronics [1]. The mass of the 
lower MACRO was ~ 4200 t, mainly in the form of boxes filled with crushed Gran Sasso 
rock. The detector had a modular structure: it was divided into six sections referred to 
as supermodules. Each active part of one supermodule had a size of 12.6 x 12 x 9.3 m 3 
and had a separate mechanical structure and electronics readout. The full detector had 
global dimensions of 76.5 x 12 x 9.3 m 3 and provided a total acceptance to an isotropic 
flux of particles of ~ 10000 m 2 sr. The total mass was ~ 5300 t. 

Data were taken from March 1989 till april 1994 with the detector under construction 
and from the middle of 1994 till the end of 2000 with the full detector. 



2 Atmospheric neutrino oscillations 

If neutrinos have non-zero masses, one has to consider the weak flavour eigenstates 
u e , v T and the mass eigenstates z/i,z/ 2 ,z/ 3 . The weak flavour eigenstates v% are linear 
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combinations of the mass eigenstates u m through the elements of the mixing matrix \J\ m : 

3 

v\ = u i™ v ™ (*) 



m=l 



In the simple case of only two flavour eigenstate neutrinos (z/ M , u r ) which oscillate with 
two mass eigenstates (1/2, z/3) one has 

f Vp = u 2 cos 9 23 + vz sin 6» 23 ^ 
\ v T = —1/2 sin # 2 3 + ^3 cos 9 2 3 

where 9 2 % is the mixing angle (9 will be used in the following). In this case one may 
easily compute the following expression for the survival probability of a beam 



-^u (t ) = l- sin 2 29 23 sin 2 [ 2 = 1 - sin 2 2fl 23 sin 2 f J (3) 

where Am 2 = m| — m 2 , and L is the distance travelled by the muon neutrino from 
production to detection. 

High energy primary cosmic rays interact in the upper atmosphere producing pions and 
Kaons, which by decay yield muons and muon neutrinos, 7r — > /jz/^; further decays of the 
muons lead to electron and muon neutrinos, /x — > v^Vil. The neutrinos are produced in a 
spherical shell at about 10- 20 km from the earth surface. Upgoing z/^'s may interact in the 
rock below MACRO or inside its lower part leading to upgoing muons, v^N — > + .... 

Upward going muons are identified by the streamer tube system (for tracking) and 
the scintillator system (for time-of-flight measurement). A rejection factor of at least 10 7 
is needed in order to separate upgoing muons from the background due to the downgo- 
ing muons. Fig. Q shows sketches of the different neutrino event topologies analyzed: 
Upthroughgoing muons, Upsemicontained (also called Internal Upgoing muons, IU), Up- 
going Stopping muons (UGS), Internal Downgoing muons (ID). The numbers of events 
measured and expected for the three topologies are given in Table H All data samples 
deviate from MC expectations; the deviations point to neutrino oscillations. 



2.1 Upthroughgoing muons 

The upthroughgoing muons come from interactions in the rock below the detector; 
the v^s have an average energy E u ~ 50 GeV [2, 3]. The data refer to a livetime of 
6.16 years (full detector equivalent). The data deviate in absolute value and in shape 
from Monte Carlo (MC) predictions, see Fig. 2. We studied a large number of possible 
effects that could affect our measurements: no significant systematic problems exist in 
the detector or in the data analyses. 

The measured data have been compared with MC simulations. For the upthroughgoing 
muon simulation, the neutrino flux computed by the Bartol group was used [4]. The 
cross sections for the neutrino interactions were calculated using the deep inelastic parton 
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Figure 2: Zenith angle distribution of upthroughgoing muons (black points). The dashed line 
is the expectation for no oscillations (with a 17 % scale uncertainty band). The solid line is the 
fit for an oscillated muon flux which yields maximum mixing and Am 2 = 2.5 • 10~ 3 eV 2 [2] [3]. 

distributions of ref. [5]. The muon propagation to the detector was done using the energy 
loss calculation in standard rock. The total systematic uncertainty on the expected muon 
flux, obtained adding in quadrature the errors from neutrino flux, cross section and muon 
propagation, is estimated to be 17 %. This uncertainty is mainly a scale error; the error 
on the shape of the angular distribution is ~ 5%. 

Fig. 121 shows the zenith angle (0) distribution of the measured flux of upthroughgo- 
ing muons. The Monte Carlo expectation for no oscillations is shown as a dashed line. 
Assuming — > v T oscillations, the best fit parameters are Am 2 = 2.5 • 10 -3 eV 2 and 
sin 2 29 = 1 with a probability of 66%; the result of the fit is the solid line in Fig. El The 
probability for no-oscillations is 0.4 %. 

Fig. shows the allowed region for — > v T oscillations in the sin 2 29 — Am 2 plane, 
computed according to ref. [6] for the upthroughgoing muon events; our region is com- 
pared with those obtained by the SuperKamiokande (SK) [7] and Soudan 2 [8] experi- 
ments. 



2.2 Matter effects, — > v T versus — ► v sterile 

Matter effects due to the difference between the weak interaction effective potential for 
muon neutrinos with respect to sterile neutrinos would produce a different total number 
and a different zenith distribution of upthroughgoing muons [9]. In Fig. EJthe measured 
ratio between the events with — 1 < cosG < — 0.7 and the events with —0.4 < cosG < is 
shown as a black point. In this ratio most of the theoretical uncertainties on neutrino flux 
and cross section cancel. The remaining theoretical error combined with the experimental 
error is estimated to be 7%. We measured 305 events with — 1 < cosQ < —0.7 and 206 
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Figure 3: 90 % c.l. allowed regions for z/„ — > z/ r oscillations from upthroughgoing muons and 
comparison with the Soudan 2 and SK allowed regions (the Soudan 2 region is now smaller). 



with -0.4 < cosQ < 0; the ratio is R = 1.48±0.13 stat ±0.10 SJ/S . For Am 2 = 2.5- 1(T 3 eV 2 
and maximal mixing, the expected value of the ratio for — > v T is R T = 1.72 while for 
Vfi — > ^ s is Rgteriie = 2.16. The maximum probabilities Pf, es t to find a value of i? r and 
of Rsteriie smaller than the expected ones are 9.4 % and 0.06 % respectively. The ratio 
of the maximum probabilities is Pbest T / Pbest sterilf , — 157, so that — > z/ s oscillations are 
disfavoured at 99% c.l. compared to the — > z/ T channel [3, 9]. 



2.3 energy estimates by MCS of upthroughgoing muons 

The oscillation probability is a function of the ratio L/E u , Eq. 3. Even if a precise 
measurement of E v is not possible, MC simulations have shown a correlation between E v 
and the residual muon energy E^. E v was thus estimated for each event by measuring 
the muon energy E^ via Multiple Coulomb Scattering (MCS) of the iduced muon in the 
absorbers. The r.m.s. of the lateral displacement for a muon crossing the apparatus on 
the vertical is omcs — 10cm/E At ( GeV). The muon energy E^ can be measured up to a 
saturation point, occurring when (Jmcs is comparable with the detector space resolution. 
A first analysis was made studying the deflection of upthroughgoing muons with the 
streamer tubes in "digital mode" . This method could reach a spatial resolution of ~ 1 cm 
which implies a maximum measurable energy of 10 GeV [10]. As the interesting energy 
region for atmospheric neutrino oscillations spans from ~ 1 GeV to tens of GeV, it is 
important to improve the spatial resolution of the detector to push the saturation point 
as high as possible. For this purpose, a second analysis was performed with the streamer 
tubes used in "drift mode" . To check the electronics and the feasibility of the analysis, 
two tests were performed at the CERN PS-T9 and SPS-X7 beams. The space resolution 
achieved is ~ 3 mm, a factor 3.5 better than in the first analysis. For each muon, seven 
MCS sensitive variables were given in input to a Neural Network (NN) previously trained 
to estimate the muon energy with MC events of known energy crossing the detector at 
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Figure 4: Ratio of events with —1 < cos9 < —0.7 to events with —0.4 < cosO < vs Am| 
for maximal mixing. The black point is the measured value, the solid line is the prediction for 
v a v r oscillations, the dash-dotted line is the prediction for v„ — * y sterile oscillations. 



different zenith angles. The method allowed to separate the upthroughgoing muons in 
4 subsamples corresponding to average neutrino energies of 12, 20, 50 and 102 GeV, 
respectively. The comparison of the 4 zenith angle distributions with the predictions of 
the no oscillations MC shows a disagreement at low energies (where there is a deficit of 
vertical events), while the agreement is restored at the higher neutrino energies. The 
distribution of the ratios R = (Data/ MC noosc ) obtained by this analysis is plotted in Fig. 
El as a function of logio(L/E u ) [3] [11]. The black points with error bars are the data; the 
vertical extent of the shaded areas represents the uncertainties on the MC predictions for 
Vy, — > u T oscillations with maximal mixing and Am 2 = 2.5 ■ 10~ 3 eV 2 . The horizontal 
dashed line is the expectation without oscillations. The black square point was obtained 
from the low energy IU sample. 

2.4 Low energy data. 

The Internal Upgoing (IU) muons come from interactions in the lower apparatus; for 
these events two scintillation counters are intercepted; the T.o.F. is applied to identify 
upward going muons. The average parent neutrino energy for these events is 4.2 GeV. For 
neutrino oscillations one expects a reduction by about a factor of two in the flux, without 
any distortion in the shape of the angular distribution. This is what is observed in Fig. 
W [12]. 

The upstopping muons (UGS) are due to external interactions yielding upgoing 
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Figure 5: Data/MC vs L/E v for upthrougoing muons (black circles) and for semicontained up-/i 
(black square). The muon energy was estimated by MCS and E u by MC methods. The shaded 
region represents the uncertainty in the MC prediction assuming sin 2 29 = 1 and Am 2 = 0.0025 
eV 2 . The horizontal dashed line at Data/MC=l is the expectation for no oscillations. 
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Figure 6: Zenith distributions for (a) the upsemicontained (IU) and (b) the upstopping plus 
the downsemicontained (UGD+ID) events. The black points are the data, the dashed lines at 
the center of the shaded regions correspond to MC predictions assuming no oscillations. The 
full line is the expectation for v„ — * v T oscillations. 
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Events 


MC-No oscillations 


R = (Data/MC noosc ) 


Up throughgoing 


809 


1122 ±191 


(0.721 ± 0.026 sto ± 0M3 sys ± 0.123^) 


Internal Up 


154 


285 ± 28 sys ± 71th 


(0.54 ± 0.04 stat ± QM sys ± 0.13th) 


Up Stop + In Down 


262 


375 ± 37 sys ± 94 th 


(0.70 ± 0.04 stat ± 0.07 sys ± 0.17th) 



Table 1: Summary of MACRO y„ — > /i events in — 1 < cos^ < (after background subtraction) 
for oscillation studies. For each topology (see Fig. ^) the number of measured events, the MC 
prediction for no-oscillations and the ratio (Data/MC-no osc) are given [2, 3, 12]. 



muons stopping in the detector. The semicontained downgoing muons (ID) are due to 
z^-induced downgoing muon tracks with vertex in the lower MACRO (Fig. 1). The two 
types of events are identified by topological criteria; the lack of time information prevents 
to distinguish the two sub-samples. An almost equal number of UGS and ID events is 
expected. In case of oscillations, the flux of the UGS muons should be reduced by 50%; 
no reduction is instead expected for the semicontained downgoing events (coming from 
neutrinos with path lengths of ~ 20 km); therefore one expects a global reduction of 25%. 

The number of events and the angular distributions are compared with MC pre- 
dictions in Table ^ and Figs. 6a,b. The data show a uniform deficit for the whole 
angular distribution with respect to predictions, ~ 50% for IU, 75% for ID + UGS; 
there is good agreement with the predictions based on neutrino oscillations with the pa- 
rameters obtained from the upthroughgoing muons. The average value of the double 
ratio R = (Data/ MC)w / (Data/ MC)id+ugs over the measured zenith angle range is 
R ~ 0.77 ± 0.07; the error includes statistical and theoretical uncertainties; R — 1 is 
expected in case of no oscillations. 

3 Search for Astrophysical Sources of High Energy 
Muon Neutrinos 

High energy u^s are expected to come from several galactic and extragalactic sources. 
Neutrino production requires astrophysical accelerators of charged particles and some kind 
of astrophysical beam dumps. 

The excellent angular resolution of our detector allowed a sensitive search for upgoing 
muons produced by neutrinos coming from celestial point sources, with a negligible atmo- 
spheric neutrino background. An excess of events was searched for around the positions 
of known sources in 3° (half width) angular bins. This value was chosen so as to take 
into account the angular smearing due to multiple muon scattering in the rock below the 
detector and by the energy-integrated angular distribution of the scattered muon, with 
respect to the neutrino direction. 

The pointing accuracy of the detector was checked with the observation of the shadow 
of primary cosmic rays (and thus on downgoing muons) by the Moon and the Sun [3] 
[13]. It was also checked with the study of the correlations between muons detected by 
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Figure 7: High energy astronomy, (a) Upgoing muon distribution in equatorial coordinates 
(1356 events), (b) The black points are the MACRO 90 % c.l. upwardgoing muon flux limits as 
a function of the declination for 42 selected sources. The solid line refers to the limits obtained 
for those cases for which the atmospheric neutrino background was zero. The limits obtained 
by the SK (open circles) and AMANDA (thin line) experiments are quoted. 



MACRO and extensive air showers detected by the EASTOP experiment located above 
the Gran Sasso massif [14]. 

In a total livetime of 6.16 y we obtained a total of 1356 events, see Fig. 7a [3, 13]. The 
90% c.l. upper limits on the muon fluxes from specific celestial sources lay in the range 
10~ 15 — 10~ 14 cm -2 s _1 , see Fig. 7b. The solid line is our sensitivity vs declination. Notice 
that we have two cases, GX339-4 (a = 255.71°, 5 = -48.79°) and Cir X-l (a = 230.17°, 
5 = —57.17°), with 7 events: in Fig. 7b they are considered as background, therefore the 
upper flux limits are higher; but they could also be indications of signals [3, 15]. 

We searched for time coincidences of our upgoing muons with sources of 7-ray bursts 
as given in the BATSE 3B and 4B catalogues, for the period from April 91 to December 
2000 [15]. No statistically significant time correlation was found. 

We have also searched for a diffuse astrophysical neutrino flux, for which we establish 
a flux upper limit at the level of 1.5 ■ 10~ 14 cm _2 s _1 [16]. 

4 Indirect Searches for WIMPs 

Weakly Interacting Massive Particles (WIMPs) could be part of the galactic dark matter; 
they could be intercepted by celestial bodies, slowed down and trapped in their centers, 
where WIMPs and anti- WIMPs could annihilate and yield — > upthroughgoing muons. 
The annihilations in these celestial bodies would yield neutrinos of GeV-TeV energy, 
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(a) (b) 

Figure 8: (a) The solid line is the MACRO upwardgoing muon flux upper limit (90% c.l.) from 
the Earth plotted vs neutralino mass m x . (b) The same as in (a) but for upwardgoing muons 
from the Sun. Each dot is obtained varying model parameters. The open circles indicate models 
excluded by direct measurements assuming a local dark matter density of 0.5GeVcm -3 . 



in small angular windows from their centers. For the Earth we have chosen a 15° cone 
around the vertical: we find 863 events. The MC expectation for atmospheric with- 
out oscillations gives a larger number of events. We set a conservative flux upper limit 
assuming that the measured number of events equals the expected ones. We obtain the 
90 % c.l. limits for the upgoing muon flux shown in Fig. |HK (it varies from 0.8 to 0.5 
10~ 14 cm -2 s _1 ) [17, 3]. If the WIMPs are identified with the smallest mass neutralino, 
our limits may be used to constrain the neutralino mass, following the model of ref. [18], 
Fig. 

A similar procedure was used to search for muon neutrinos from the Sun, using 10 
search cones from 3° to 30°. In the absence of statistically significant excesses the muon 
flux upper limits are at the level of 1.5 — 2 • 10~ 14 cm -2 s . The limits are shown in Fig. 
o D as a function of the WIMP (neutralino) mass. 



5 Neutrinos from Stellar Gravitational Collapses 

A stellar gravitational collapse (GC) of the core of a massive star is expected to produce 
a large burst of all types of neutrinos and antineutrinos with energies of 7 — 30 MeV and 
with a duration of < 10 s. The zVs can be detected via the process v e + p — > n + e + in the 
liquid scintillator. About 100 -j- 150 v e events should be detected in our 580 t scintillator 
for a stellar collapse at the center of our Galaxy. 

We used two electronic systems to search for bursts of P e 's from stellar gravitational 
collapses. The first system was based on the dedicated PHRASE trigger, the second one 
was based on the ERP general trigger. Both had energy thresholds of ~ 7 MeV and 
recorded pulse shape, charge and timing informations. Immediately after a trigger, the 
PHRASE system lowered its threshold to 1 MeV, for a duration of 800/i s, in order to detect 
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(with a ~ 25 % efficiency) the 2.2 MeV 7 released in the reaction n + p — > d + 72.2MeV 
induced by the neutron produced in the primary process. 

A redundant supernovae alarm system was in operation, alerting immediately the 
physicists on shift. We defined a general procedure to alert the physics and astrophysics 
communities in case of an interesting alarm. Finally, a procedure to link the various 
supernovae observatories around the world was set up [3, 19]. 

The effective MACRO active mass was ~ 580 t; the live-time fraction was ~ 97.5 %; 
we kept always ruming at least 1-2 supermodules. No stellar gravitational collapses were 
observed in our Galaxy from the beginning of 1989 to the end of 2000 [3] [19]. 

6 Conclusions 

The MACRO detector took data from 1989 to the end of year 2000 and it obtained 
important results in all the items listed in the proposal, in particular: 

• Atmospheric neutrino oscillations. Analyses of different event topologies at different 
energies, the exploitation of the muon Coulomb multiple scattering in the detector 
give strong support to the hypothesis of — >• v T oscillations with Am 2 = 0.0025e^ 2 
and maximal mixing. The statistical and systematic significance of our global data 
is well above at the 5 standard deviation level. We are at present performing some 
global analyses using also the new data on atmospheric muon flux measurements. 

• High energy muon neutrino astronomy. MACRO has been highly competitive with 
other underground experiments thanks to its large acceptance and good angular 
accuracy. It has been limited by its livetime and the size of the detector. 

• Search for bursts of v e from stellar gravitational collapses. MACRO was sensitive 
to supernovae events in our Galaxy, it started the SN WATCH system, and for a 
certain time it was the only detector in operation. 

• Sensitive indirect searches have been carried out for possible Dark Matter candidates 
like the neutralinos, looking for upgoing muons from the center of the Earth and 
from the Sun. 

It may be pointed out that MACRO also obtained the best flux upper limits for GUT 
Magnetic Monopoles over the widest j3 range [20]. 
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